Breakdown of adiabatic transfer schemes in the presence of decay 
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In atomic physics, adiabatic evolution is often used to achieve a robust and efficient population 
transfer. Many adiabatic schemes have also been implemented in optical waveguide structures. 
Recently there has been increasing interests in the influence of decay and absorption, and their 
engineering applications. Here it is shown that contrary to what is often assumed, even a small 
decay can significantly influence the dynamical behaviour of a system, above and beyond a mere 
change of the overall norm. In particular, a small decay can lead to a breakdown of adiabatic transfer 
schemes, even when both the spectrum and the eigenfunctions are only sightly modified. This is 
demonstrated for the decaying version of a STIRAP scheme that has recently been implemented in 
optical waveguide structures. It is found that the transfer property of the scheme breaks down at a 
sharp threshold, which can be estimated by simple analytical arguments. 
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Introduction. The adiabatic theorem of Hermitian quan- 
tum mechanics provides the basis for many experimen- 
tal schemes to manipulate and control quantum systems. 
Prominent examples from atomic physics include the 
RAP and STIRAP schemes [1 for population transfer 
in effective two or three level systems. In many applica- 
tions, however, quantum systems are open, that is, their 
states have finite lifetimes described by complex energies. 
These more realistic situations cannot be described by 
Hermitian quantum mechanics, but require an additional 
anti- Hermitian part in the Hamiltonian (see, e.g., ^ i3j 
and references therein). Due to the analogy of the time 
dependent Schrodinger equation and the paraxial approx- 
imation for light propagation in optical media [H |5] , the 
dynamics generated by non-Hermitian Hamiltonians can 
conveniently be realised using optical systems [6]. 

It has previously been pointed out in the literature that 
the adiabatic theorem does not necessarily hold for non- 
Hermitian systems 0, and it has recently been shown 
that this can lead to new effects in the "deep" non- 
Hermitian regime in the presence of exceptional points 
[8 . However, it is often assumed that small decay rates 
do not modify the systems behaviour drastically, and 
their effects on adiabatic behaviour seem hitherto not 
to have been fully appreciated. 

Here we show that even a small decay rate, which does 
not modify the static behaviour of a system significantly, 
can lead to a breakdown of adiabatic transfer properties. 
We demonstrate this fact for a STIRAP-related scheme, 
which is readily implemented in optical waveguide struc- 
tures. In this scheme a sharp threshold for the breakdown 
of the transfer property is observed, which we calculate to 
a good approximation using simple analytical arguments. 
Breakdown of adiabaticity in non-Hermitian systems. 
The adiabatic theorem states that a system initially pre- 
pared in an eigenstate remains in the corresponding in- 
stantaneous eigenstate if the system parameters are var- 



ied infinitely slowly and the corresponding energy level 
is nondegenerate at all times. For finite parameter varia- 
tions, there are small transition amplitudes between the 
adiabatic states, which are typically of order 0(e), where 
e denotes the small adiabatic parameter [9l [10] . In the 
case of complex eigenvalues, however, the amplitudes of 
the instantaneous eigenstates are themselves exponen- 
tially decreasing in time. This can lead to a situation in 
which the small nonadiabatic transition amplitude grows 
exponentially in time relative to the adiabatic amplitude, 
if the adiabatic state is not the one with the smallest de- 
cay rate. In other words, the effect is caused by the 
dominance of the single gain mode of the time evolution 
operator [8\ We will now demonstrate how this general 
phenomenon can lead to a breakdown of the STIRAP 
scheme when the final state has a finite lifetime. 
The STIRAP scheme. We begin with a brief summary of 
the conventional STIRAP scheme for a three-level quan- 
tum system described by a time dependent Hamiltonian 
of the form [H [11] 

/ v{t) \ 
H{t) = v(t) w(t) . (1) 
\ w(t) ) 

The instantaneous energies of the system are given by 

^0 = 0, = (2) 

with the corresponding eigenstates 

(cos6> \ -| / sin^ \ 

, 1^^) = -^ ±1 , (3) 
- sin i9 ) v2 y cos6> / 

where tan6> — vjw. 

A variation of the coupling parameters v{z) and w(z) 
fulfilling ''^^l^^^'-'X oo and ""^Z/^""^^, rotates the 
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FIG. 1. The left figure shows the eigenvalues of the Hamilto- 
nian ([T]), with the parameter dependence Q for a = 5. The 
right figure shows the squared absolute values of the compo- 
nents of the adiabatic eigenstate \(p^) (first, second, and third 
component shown in blue, green, and red, respectively). 



instantaneous eigenstate \(p^) from |3) to |1). Thus, adia- 
batic parameter variation leads to a complete population 
transfer from |3) to |1). Furthermore, since the second 
component of \(p^) is identically zero at all times, the 
level 1 2) is never significantly populated during the pro- 
cess. The famous STIRAP scheme of atomic physics is 
an example of this adiabatic population transfer scheme 
[T]. In this context, while the states |1) and |3) are usu- 
ally assumed stable, the intermediate state |2) generally 
couples to a continuum of other states, leading to a de- 
cay. Due to the vanishing population in |2), however, the 
scheme is insensitive to this decay [1 . Surprisingly, what 
appears not to have been noticed so far is that a decay 
from one of the other states can lead to a breakdown of 
the scheme. In what follows we shall demonstrate this for 
a decay from the final state, which is in practice often an 
excited energy state and thus less stable than the initial 
state. 

Recently, STIRAP-type schemes have also been imple- 
mented in optical waveguides [12^ , where the prop- 
agation distance z takes the role of time. The optical 
realisation allows for a straightforward experimental im- 
plementation of an additional decay of varying strength 
using absorbing materials. Thus, in what follows we shall 
proceed our analysis in the waveguide context. In partic- 
ular, we study a generalisation of the setup investigated 
in p2J, which consists of two parallel (left and right) 
waveguides with an additional diagonally directed cen- 
tral waveguide. Since the couplings between neighbour- 
ing waveguides depend exponentially on their distances, 
this system can be described by a Hamiltonian of the 
type ([T]) with couplings of the form 



v{z) = l/^(z) =e-^("-^/2)/^. 



(4) 



The instantaneous eigenvalues of this system are depicted 
in the left panel in Fig. [l] The right panel shows the com- 
ponents of the instantaneous eigenstate corresponding to 
the zero energy eigenvalue. For sufficiently large values 
of a and L, a complete intensity transfer between the 
right and the left waveguide is achieved without ever sig- 
nificantly populating the central waveguide, as has been 
experimentally demonstrated in [12j. 
STIRAP with decay in the final state. Let us now in- 




FIG. 2. Real (top) and imaginary (middle) parts of the en- 
ergies for 7 = 0.1 (left) and 7 = 0.5 (right), and a = 5. 
The solid black line shows the numerically obtained eigenval- 
ues, and the dashed red line those obtained from first order 
perturbation theory. The solid lines in the lower panel show 
the components of the adiabatic eigenstate Icp'^) (colours as in 
figure 1), the dashed lines show the perturbative result. 



vestigate how absorption in the target (left) waveguide 
modifies this transfer scheme. The additional decay is 
modelled by an imaginary energy —ij in the first diago- 
nal element of the Hamiltonian ([T]). The eigenvalues for 
small 7 can be obtained via first order perturbation the- 
ory that leaves the real parts unaltered and yields the 
additional imaginary parts: 



lm{Eo) ^ -7cos^(9 
lm{E±) ^ ^^^^ ^ 



-77 



2(t;2 + u,2)- 



(5) 



The eigenstate corresponding to Eq^ in first order pertur- 
bation theory, is given by 



cos 

\(p^) = I ij cos 6 sin / ^/v^ 
— sin^ 



(6) 



In the upper and middle panels of Fig. |2] we show the 
numerically obtained energy levels in comparison to the 
first order perturbation theory for a = 5, and two dif- 
ferent values of 7 as a function of z/L. The energies are 
well described by the perturbative equation. In the lower 
panels of the figure the components of the instantaneous 
eigenstate \(f^) are shown. As in the Hermitian case, this 
state populates the right waveguide for z = 0, while it 
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FIG. 3. Numerically obtained transfer probability (P = 1 
is red, P = is blue) as a function of L and 7 for a = 5 
(left) and a = 8 (right) . The dashed white line shows the 
analytical transition boundary based on the amplification of 
the nonadiabatic coupling for the Hermitian system using the 
approximation (17), the solid black line is the semianalytical 
critical boundary using the numerically obtained nonadiabatic 
coupling for the non- Hermitian system with 7 = 0.2. 



is localised in the left ioi z = L. However, the central 
waveguide is also partly populated for intermediate val- 
ues of z. 

Assuming that the adiabatic theorem holds, one would 
conclude from the behaviour of the eigenvalues and eigen- 
states that the STIRAP scheme should lead to a full 
transfer of the population (up to an overall decay) for 
moderate values of 7. However, a numerical simulation 
of the total transferred probability, as a function of 7 
and L, yields a different result, as depicted in Fig. |3] 
There is a sharp transition for relatively small values of 
7 at which the STIRAP scheme breaks down. This is 
related to the relative growth of the nonadiabatic transi- 
tions mentioned above. In what follows we shall discuss 
in detail how this transition boundary emerges. 
Adiabatic transition boundary. The transition probability 
to the left waveguide is defined by 



P = 



\{mm 



(7) 



Here we consider the transition probability starting ini- 
tially in the state |(/:?^), corresponding approximately to 
the right waveguide, that is, 



1^(0)) = 1/(0)) -|3). 



(8) 



Let us express the solution of the time dependent 
Schrodinger equation as a sum of adiabatic and non- 
adiabatic parts in the instantaneous eigenbasis: 

\^{Z)) = ^ad{z)\^'{z)) + Yl ^nonadizWiz))^ (9) 

where the nonadiabatic coefficient of the instantaneous 
states \(p^) are equal due to the symmetry of the problem. 
On account of the relation \(p^{L)) ^ |1), the transfer 
probability in Eq. ^ can be estimated as 



P = 



\l^ad{L)\^^2\l^nonad{LW 



(10) 



In the Hermitian case the transfer is successful as long 
as the nonadiabatic term is negligible, which is the case 
for sufficiently large values of L. However, in the non- 
Hermit ian case, we have to take into account that the 
adiabatic states have a finite lifetime, that is, they decay. 

To a good approximation the transition between the 
instantaneous states occurs in the vicinity of the avoided 
crossing of the energies, close to z = L/2 (see Fig. [2|. 
The states \(p^) decay only slowly after this point, since 
the imaginary parts of their energies are small. Thus we 
find 



1^. 



'nonad 



(L)| « VP, 



nonadi 



(11) 



where Pnonad denotes the nonadiabatic transition prob- 
ability due to the finite parameter variation. The state 
|(/?^), which we wish to follow adiabatically, on the other 
hand, has a small imaginary energy and is approximately 
stable before the nonadiabatic transition takes place, 
while the remaining population decays after the transi- 
tion. Thus we estimate 

\ll^ad{L)\ = \/l - 2Pnonad exp ^ j IvCiE^dz^ 

- v/l-2Pnonade-^^/^ (12) 

where the integral is calculated using the eigenvalues in 
first order perturbation theory in equation The fac- 
tor of 2 in front of the nonadiabatic transition probability 
accounts for transitions into the two states \^^). 
The population transfer is successful if 



\^ad{L)\^:>'^\^nonad{L)\^- 



(13) 



Treating ( |13[ ) as an equality and using ([11} and ([12|, we 
obtain the threshold value 



In 



1 



2P 

All ri. 



1 /L. 



(14) 



In the same way, the characteristic width of the threshold 
can be estimated by 5^cr ~ 2/L. Since 5^cr ^ "fcr for 
small Pnonad Eq. (14), the observed threshold appears 
to be sharp. 

Since we expect 7cr to be relatively small to get an esti- 
mate, we approximate the nonadiabatic transition prob- 
ability with its value in the Hermitian case where 7 = 0. 
Unfortunately, even in this case no analytic expression 
for Pnonad IS kuowu [14] . Usiug a Landau-Zener type 
approximation, however, it can be estimated as 



nonad 



exp ^-2Im j \e^- EQ)d^ , (15) 



where zq denotes the position of the exceptional point 
nearest to the real z-axis [15 . 

In the Hermitian case exceptional points of all three 
levels (EP3) appear for Ei = E2 = E^, = correspond- 
ing to complex values of z given by the equation v = :^iw. 
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FIG. 4. Nonadiabatic transition probabilities as a function 
of L for a = 5 for different values of 7, 7 = (light blue), 
7 = 0.1 (magenta) and 7 = 0.25 (black). For comparison also 
the estimation based on equation ( |17| ) is shown (black dashed 
line) . 



Using Q, we find that they are located at 

= Q + ^(1 + 2n)^^ , ne Z. (16) 

The EP3 with the smallest positive imaginary part of z 
is zq. From Eq = and E^{L/2 + i^) = Vv^Tw^ 



Y^2 cos(2a^/L), we can compute the integral in (15) as 

plm Zq 



exp 



exp 



V'2 cos{2a^/L)d^ 



(17) 



Substituting this expression into (14) yields an analytic 



estimate for the critical decay rate, 



I cr 



L 



(18) 



which for large values of L simphfies to 



..LZ 



ln(2) 



(19) 



In Fig. |3] this value is shown as a white dashed line. 
It can be seen that it is a good estimate for the exact 
transition boundary. It fails, however, to describe a slow 
decrease of jcr for large values of L. This decrease is due 
to the fact that the Landau- Zener approximation for the 



nonadiabatic corrections assumes a vanishing coupling at 
the beginning and the end of the time evolution, which 
is not the case in the STIRAP scheme investigated here. 

The exact nonadiabatic transition probability can be 
obtained by a numerical integration. We show the re- 
sult for three different values of 7 in a semilogarithmic 
plot as a function of L in Fig. [4j The light blue line 
corresponds to the Hermitian case. For comparison, the 
Landau-Zener estimate (17) is also shown (dashed black 



line). The initially approximately exponential decrease 
of the transition probability with increasing L is well de- 
scribed by the Landau-Zener type expression. At a crit- 
ical value of L, however, the transition probability be- 
gins to oscillate around an only slowly decreasing mean 
value. Similar behaviour is typically observed in STIRAP 
schemes [14] . The presence of decay leads to a smoothing 
of the oscillations, and an increase of the critical value of 
L. It can also be seen that the slope in the region of ex- 
ponential decay is slightly decreased by the decay. This 
explains why the white lines in Fig. [3] slightly overesti- 
mate the critical value of 7. As expected the quality of 
the approximation is better for the case a = 8, where the 
boundary is located at smaller values of 7. 

We can obtain an excellent approximation for the 
transfer boundary by using the numerical values of the 
nonadiabatic transition probability for a fixed value of 
7 = 0.2, close to the actual boundary, as an input in 
equation (14). This is demonstrated in Fig. |4j where the 
solid black lines depict the thus obtained result. 
Conclusion. We have demonstrated that even a tiny 
decay rate can significantly influence the dynamical 
behaviour of a system, in particular in the context of 
adiabatic time evolutions. This is due to a competition 
between small nonadiabatic transition amplitudes and 
relative exponential growths of the decaying adiabatic 
eigenstates. In particular, we have shown for a STIRAP 
related scheme, which can be implemented straight- 
forwardly using optical waveguides, that the adiabatic 
transfer behaviour breaks down at a sharp threshold 
for relatively small decay rates. The critical value of 
the decay rate has been estimated by simple analytical 
arguments. 
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